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ABSTRACT 

We present a detailed study of the star cluster population detected in the galaxy NGC 922, one of 
the closest collisional ring galaxies known to date, using HST/WFPC2 UBVI photometry, population 
synthesis models, and N-body/SPH simulations. We find that 69% of the clusters are younger than 
7Myr, and that most of them are located in the ring or along the bar, consistent with the strong Ha 
emission. The cluster luminosity function slope of 2.1-2.3 for NGC 922 is in agreement with those 
of young clusters in nearby galaxies. Models of the cluster age distribution match the observations 
best when cluster disruption is considered. We also find clusters with ages (>50Myr) and masses 
(>1O 5 M ) that are excellent progenitors for faint fuzzy clusters. The images also show a tidal plume 
pointing toward the companion. Its stellar age from our analysis is consistent with pre-existing stars 
that were stripped off during the passage of the companion. Finally, a comparison of the star-forming 
complexes observed in NGC 922 with those of a distant ring galaxy from the GOODS field indicates 
very similar masses and sizes, suggesting similar origins. 

Subject headings: galaxies: individual (NGC 922) — galaxies: star clusters — galaxies: interactions. 



1. INTRODUCTION 

Collisional ring galaxies represent a very specific case of 
galaxy interaction for which a compan ion drops through 
a mu c h more massive sp i ral ga l axy (e.g.lLvnds &: Toomrel 
[19761: iThevs fc Spiegel II9771 IHernauist fc Weill I1993H . 
First it requires very specific initial conditions to occur. 
The mass of the companion must be small compared to 
the main galaxy and the primary must be a disk galaxy. 
Second, the characteristic ring of star formation only 
lasts for a few 10 s yr, which is not a large temporal win- 
dow compared to the Hubble timescale. Although they 
are rare objects in the local univers e, they are thought 
to be more common at high redshifts ()Laverv et al.lll996t 
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Elmegreen & Elmegreen 2006). 

NGC 922 is a drop-through ring galaxy for which the 
perturber, respo nsible for the ring pattern, has recently 
been discovered (jWong et al.l 120061 hereafter W06). Us- 
ing simulations, the authors found that NGC922's 
companion (named S2) interacted gravitationally with 
NGC 922 about 330 Myr ago through an off-centered 
collision, creating the C-shape morphology seen in the 
disk of NGC 922. The dwarf compact companion is 
about 20 times less massive than NGC 922, the latter 
having an estimated mass of 2.8xl0 10 M Q . NGC 922 
is one of the nearest colli sional ring galaxies known, 
at a distance of 43Mpc (Me urer et al.l 120061) . about 
thre e times closer than the iconic Cartwheel Galaxy 
(e.g.lFosbury fc HawardenH l977t ID avies fc Morton 1982; 



(e.g. .tosburv &; riavv 
IStruck et alT li996 N ). 
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It has a heliocentric velocity of 
3077 km s" 1 , and an estimated star formation rate (SFR) 
of 7-8 M Q yr -1 . A more detailed description of the phys- 
ical properties of both NGC 922 and S2 can be found in 
the work of W06. 
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Fig. 1— Left: Three-color HST/WFPC2 image of NGC922. Blue: F300W; Green: F547M; Red: F814W. The image is 2.5'x2.0' 
(~30.7x24.6 kpc 2 ). Right: Ha image of NGC922 from SINGG. The resolution is ~1.5", lower than the HST image. Blue: R-band; Green: 
Ha+ continuum; Red: Ha with continuum substracted. North is up and east is left. 



NGC 922 is a very interesting object to study. Since the 
companion responsible for the star formation activity is 
known, dynamical simulations can be better constrained 
and its recent star formation history can be modeled. 
Furthermore, because of its proximity, its stellar clus- 
ter population can be resolved and observed with the 
Hubble Space Telescope (HST). In this paper, we present 
recent HST optical images of NGC 922. The overall goal 
is to study the recent star cluster formation history of 
NGC 922 and to compare it to the dynamical models in 
order to constrain the models and the recent dynamical 
history of NGC 922 as well as better understanding the 
star formation and cluster evolution processes. 

In the next section, we present the HST data as well 
as the photometric analysis of the star clusters detected 
in NGC 922. In section [3] we present the synthesis mod- 
eling performed to estimate cluster ages and masses. In 
section |4] we describe the star cluster properties. Sec- 
tion [5] discusses the effect of cluster dissolution on the 
observed star cluster population. We then present an N- 
body/SPH simulation of the collision and compare the 
results to the observations (i}6]). In section [7] we investi- 
gate the large-scale star formation of NGC 922 and com- 
pare it to distant ring galaxies. Section [8] presents diffuse 
features detected in NGC 922. We then discuss the fate 
of NGC 922 and its star clusters in section HI We con- 
clude in gTS 

Six sets of images of NGC 922, two slightly different 
pointings with three exposures each, have been obtained 
on 25 and 26 June 2007 using the Wide Field Plane- 
tary Camera 2 (WFPC2) onboard HST in the four filters 
F300W, F439W, F547M, and F814W (hereafter U, B, V, 
and I for simplicity). The F439W and F547M filters were 
chosen over F450W and F555W to avoid contamination 
by the nebular emission lines of [O in] and H/3 in the pass- 
bands. The WFPC2 field covers the whole ring galaxy, 
but excludes S2. 

The data were first reduced by the STScI calibration 
pipeline, which includes treatment for the bias, dark, and 
flat fields. Then we corrected for warm pixels, geometric 



distortion and flux calibration, as well as charge transfer 
efficiency. We used the variable-pixel linear reconstruc- 
tion method (or drizzling) and the IRAF task multidriz- 
zle to combine the individual exposures, remove cosmic 
rays and to improve the spatial resolution of the final im- 
ages. The chips were drizzled individually to an output 
pixel size of 0.075" pixel" 1 (15.6 pc pixel -1 ) for the WF 
chip. 

A three-color image from our HST/WFPC2 data is 
shown in Figure [T] together with a Ha image from the 
Survey for Ionizati on in Neutral-Gas Galaxies (SINGG; 
IMeurer et alj |2"006) . Figure [5] presents the main morpho- 
logical features of NGC 922 over-plotted on the F814W 
WFPC2 image. The blue ring has a projected diameter 
of 15 kpc and is particularly prominent on the eastern 
side, giving a "C" like morphology to NGC 922. While 
the eastern side of the ring displays sharp edges with 
star-forming regions well distributed along it, the west- 
ern side shows a patchy distribution and weak diffuse Ha 
emission. A 9 kpc off-center bar, which appears yellow 
in Fig. QJi, displays spokes emanating from it. A rel- 
atively isolated star-forming region is also observed on 
the South- West portion of the ring but is disconnected 
from the "C" . Another star-forming structure, like an 
accent over the "C" , is observed north of the ring, and 
merges quite well with the ring near the bar end. On the 
Ha image an inner arc can be seen internal to the ring, 
perhaps as an extension of the "accent". However it is 
impossible to conclude here if this structure is physically 
real or just a random distribution of star-forming regions 
misleading the eye. 

To study the star cluster properties, we performed 
PSF photometry using the IRAF/DAOPHOT package. 
PSF photometry is required because of crowding. For 
each chip and each pointing, we used bright unsaturated 
sources to create a PSF model of 4.5 pixels radius and 
a FWHM~2.5 pixels, which slightly varies depending on 
the chip. In the case of the PCI chip, there were too few 
sources to create a reliable PSF model. For this reason, 
and because the same region was observed on the WF2 
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Fig. 2. — Morphological features of NGC 922 over-plotted on the F814W image. The circle shows the expanding ring of star formation. 
The + symbol identifies the nucleus. The white lines follow the bar. The dotted lines indicate the location of possible spokes. The polygon 
highlights the tidal plume. The solid line, North of the circle, shows a feature of star formation outside of the expanding ring. The physical 
scale is the same as in Fig. [l] 



chip of the other pointing, we are using only the photom- 
etry based on the WF chips for this work. A summary of 
the observational datasets is presented in Table [T] The 
absolute photometric calibration is in the Vegamag sys- 
tem and was obtained with the zero points calculated for 
each chip as described on the STScI/WFPC2 cookbook 
websit^E Because the Multidrizzlc software converts the 
WFPC2 images from DN to e- in order to properly re- 
ject cosmic rays, we adjusted the photometric zeropoint 
to account for this gain. We calculated and applied the 
aperture correction for each of the four filters and chips. 

To determine the photometric uncertainties and the 
completeness level of each cluster, we generated a list of 
artificial sources of various magnitudes and background 
levels for each combination of chip and filter. This is 
critical to take into account that more clusters will re- 
main undetected in regions of high level background (of- 
ten linked to crowding), and that the photometry will 
also be less accurate. We used this list as input for the 
DAOPHOT/addstar task to generate aritificial sources 
in the WFPC2 images. We then performed the PSF 
photometry again on the new images with added stars. 
We used the results to create functions of photometric 
uncertainties and completeness levels that takes into ac- 
count the background level. Then we used these func- 
tions to associate the proper photometric uncertainty 
and completeness level values to each detected source 
and each filter. On average, the 50% completeness lim- 
its are reached at 23.8, 24.3, 26.2, and 27.1 magnitudes 
for the U, B, V, and I filters respectively. We also ap- 
plied a foreground ex tinction correction of E(B— V)=0.02 
(jSchlegel etaLlll998lh We find that these detection lim- 



its may include individual bright stars among faint clus- 
ter samples (see discussion in H4.2[) . Finally we do not 
expect much contamination from background galaxies, 
i.e. about 5 to 15 galaxies over the entire field, to 
our exposure depth , based on the Hubble Deep Field 
(|Shanks et al.llT998h . 

We created four lists for all the clusters detected in 
each of the four filters. We then cross-correlated the 
four lists to find sources with matching coordinates. The 
matching objects had to be within 2.5 pixels from each 
other to be considered the same physical source. Overall 
we observed 2248 star cluster^ detected in the four fil- 
ters simultaneously through PSF photometry (hereafter 
UBVI sample). An additional 1569 clusters were also 
observed in the three BVI filters (hereafter BVI sam- 
ple), but displayed no conterpart in the U-band image. 
Because of the given spatial resolution of 39 pc (2.5 pix- 
els), it is probable that some of the sources detected con- 
tain multiple clusters. A completeness correction was 
applied to each cluster found in each filter and based on 
the completeness level determined in that filter. For si- 
multaneous detections in multiple filters, we combined 
the completeness corrections of the appropriate filters. 

An examination of the residual images showed that 
three extended sources were significantly bigger than the 
PSF size, and were therefore badly fitted. They comprise 
two large clusters and the nucleus. For those sources, we 
performed a simple aperture photometry using polygons 
and the IRAF task polyphot. The photometric results for 
those three sources are given in Table [5] 

3. STELLAR POPULATION MODELING 



3 http:/ /www. stsci.edu/hst/wfpc2/analysis/wfpc2_cookbook. html 



4 Due to the very large number of clusters, the detailed photom- 
etry is not published here, but is available upon request. 
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TABLE 1 

HST/WFPC2 Observational Data Parameters of NGC922 



Filter 


Pointing 


Exp. Time 


Detection Limit 


Chip 


PSF FWHM 






sec 


mag 




pixels 


F300W 


A 


1100 


23.8 


WF2 


2.5 




A 






WF3 


3.3 




A 






WF4 


2.6 




B 


1100 




WF2 


2.5 




B 






WF3 


2.8 


F439W 


A 


600 


24.3 


WF2 


2.3 




A 






WF3 


2.5 




A 






WF4 


2.3 




B 


600 




WF2 


2.3 




B 






WF3 


2.5 


F547M 


A 


600 


26.2 


WF2 


2.3 




A 






WF3 


2.5 




A 






WF4 


2.5 




B 


600 




WF2 


2.1 




B 






WF3 


2.1 


F814W 


A 


600 


27.1 


WF2 


2.2 




A 






WF3 


2.8 




A 






WF4 


2.3 




B 


600 




WF2 


2.5 




B 






WF3 


2.5 



We used the populat ion synthesis model Starburst9£0 
(|Leitherer et al.l 119990 to estimate cluster ages and 
masses. We used an instantaneous burst of sta r forma- 
tion w ith a standard initial mass function (IMF; iKroupal 
2002) ranging from 0.1 to 100 M Q , as well as the Padova 
stellar evolutionary tracks with stars from the asymptotic 
giant branch. A solar metallicity was chosen to match the 
result of W06 and based on the nuclear spectrum from 
the 6dF Galaxy Redshift Survey (pones et al.ll2004[ ). For 
comparison, we also performed the synthesis using dif- 
ferent model properties with 1) an IMF ranging from 0.1 
to 30 M Q and 2) with a metallicity of 0.4 Z . However, 
considering the photometric uncertainties, these two ad- 
ditional models did not give significantly different results. 

For an accurate comparison between the stellar popu- 
lation models and our photometric data, we applied the 
IRAF/SYNPHOT package on the Starburst99 spectra to 
calculate the theoretical magnitudes that would be seen 
in our four WFPC2 filters. Then we performed a \ 2 fit- 
ting on the photometric data points, basically the slope 
of the spectral energy distribution (SED), to obtain the 
best-fit age. Finally, we calculated the best-fit mass for 
each cluster by averaging the magnitude differences be- 
tween the observation and the models in the four filters 
for the UBVI sample, and in the three filters for the BVI 
sample. The mass and age uncertainties were estimated 
by fitting the models within the photometric error bars 
using the bootstrap technique. 

Most of the clusters found in this galaxy are very 
young, with a mean age of 16 Myr and a median age 
of 6 Myr for the UBVI sample, and a mean age of 34 Myr 
and a median age of 6 Myr for the BVI sample. Inter- 
estingly, the sources of the BVI sample are mainly faint 
sources and we obtained young ages for most of them. 

3.1. Photometric Modeling and Extinction 

In an effort to estimate the extinction for each clus- 
ter, we used the Q- method applied to cluste r col- 
ors as described by Uohnson fc Morgan! ([1953) and 

5 http:/ /www. stsci.edu/science/starburst99/ 



Ivan den Bergh fc Hagenl (|1968t ). This method allows, in 
principal, an accurate age estimate of the clusters since 
the Q-Parameters are independent of reddening (see a 
good example of this method applied to star clusters by 
iWhitmore et al.l I1999T ) . However, this method requires 
very good photometry because the error bars get larger 
when color information is involved. For this reason, and 
despite the high quality of the WFPC2 images, we did 
not obtain extinction estimates that are precise enough 
to apply to each cluster. The ages obtained with the 
Q-Parameter method were overestimated by ^100 Myr 
compared to the more precise SED fitting described pre- 
viously, which biased the extinction calculation. Note 
that the overall E(B-V) values we obtained were consis- 
tent with the dust features seen by eye in the WFPC2 
images, with more obscured zones where larger E(B-V) 
values were found, but the cluster-to-cluster extinction 
variation was too large to reliably apply it to each clus- 
ter. Although most clusters found in this work are al- 
ready very young and cannot be that much younger, we 
must keep in mind that we probably overestimated the 
age of some clusters, especially the reddest ones. 

4. THE STELLAR CLUSTER POPULATION 

4.1. Cluster Luminosity Function 

Figure |31 shows the cluster luminosity functions (CLF). 
For each filter (one per panel) we compare the CLF of 
four groups. The first group contains all the sources de- 
tected in each individual filter and is represented by the 
black line. The second group, in blue, contains only the 
clusters from the UBVI sample. We can see that many 
faint sources are lost when applying the four color detec- 
tion criteria. The third group, in red, contains the clus- 
ters from the BVI sample, i.e., those undetected in the 
U-band. Several reasons can explain the non-detection 
of a cluster in the U-band: a) the cluster is old, b) the 
cluster is faint or c) the cluster is reddened by dust. The 
CLF in Figure |3] clearly shows that most of the sources 
undetected in the U-band are simply faint, and there- 
fore, probably simply less massive. This is consistent 
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TABLE 2 

Large Sources with Aperture Photometry 



Source 


RA 

degrees 


DEC 

degrees 


F300W 
mag 


F439W 
mag 


F547M 
mag 


F814W 
mag 


Age 
Myr 


Log(M steiiQr ) 
M 


Area 
kpc 2 


Cluster #1 
Cluster #2 
Nucleus 
Tidal Plume 


36.27565 
36.27287 
36.26851 


-24.79476 
-24.79845 
-24.78820 


18.4±0.4 
17.2±0.3 
15.34±0.04 


18.6±0.4 
17.4±0.3 
14.88±0.01 
17.0 


19.8±0.4 
18.8±0.3 
16.08±0.02 
15.06 


19.7±0.5 
18.5±0.4 
15.56±0.03 
15.45 


g+0 
°-2 
g+0 
°-2 

qn+° 
yu -20 

1300+ 3700 
iouu_ 300 


5.5±0.2 

6 0+ ' 1 

D.U_ Q 2 
« 1+° 

8.9±0.2 


41.1 
143.9 
328.3 
30260.3 



Note. — The RA and DEC positions arc in J2000 and approximations and are not necessarily related to the exact photometric center of 
the sources. 



with the young mean and median ages found for this 
sample as noted in section [3] This is also confirmed by 
the fact that 59% of the clusters from the BVI sample 
have masses smaller than 7x10 3 Mq. The fourth curve 
in cyan is the CLF for the UBVI and BVI samples com- 
bined together. We can see that this combined sample 
is in better agreement with the raw sample (black) at 
fainter luminosity compared to the UBVI sample alone 
(blue). 

We calculated the sl o pe a of the CLF as described by 
lElmegreen fc Efremovi (|1997h . For the raw samples, we 
obtained slope values of 2.15±0.08, 2.22±0.05, 2.05±0.05 
and 2.10±0.05 for the U, B, V, and I bands respec- 
tively. The slopes tend to steepen to ~2.3 for clusters 
with UBVI simultaneous detection, when we fit up to 
the 50% completeness limit magnitude. Our values are 
consistent wi th those of young clusters found in nearby 
galaxies (e.g. iWhitmore et alJll999t lLarse 3 120021 ,. 

4.2. Age Distributions 

In Figure 21 we present spatial maps of clusters for 
various age ranges. The age group 0-7 Myr corresponds 
to clusters capable of ionizing their surrounding gas and 
therefore producing Ha emission. They represent 69% 
of the detected clusters, in number, and they are mostly 
located within the ring of star formation and along the 
bar, as well as in a more isolated region, South- West of 
the nucleus. Their location is consistent with the Ha 
emission (Fig. [TJ. We also find that 19% of the clusters 
have ages between 7 and 50 Myr, and only 3% are older 
than 100 Myr. The oldest clusters are mainly located in 
the nuclear region, but are also seen along the bar and 
in the disk. One must keep in mind that some of them 
might be reddened and therefore younger than estimated 
here. 

Figure [5] shows the magnitude-age distributions of 
star clusters in NGC 922. It can be seen that at 
older ages (> 100 Myr), only the most massive clusters 
(>1x10 5 Mq) can be detected in this galaxy. They are 
not numerous and they mostly have just BVI photome- 
try. Those sources must be analyzed with care since they 
are mostly related to the nuclear region, or other very 
dense regions where the field is very crowded and lumpy, 
and therefore may include underlying stellar populations. 
Also, some of them may be particularly affected by ex- 
tinction and may have ages biased toward larger values, 
which would lead to overestimated masses. 

W06 estimate the impact epoch to be ^300 Myr ago. 
Because of the small number of clusters at similar ages 
and their large age uncertainties, we do not believe it 



is possible to constrain the impact epoch from the clus- 
ter age distribution with the current dataset. Around 
10-30 Myr, there is an effect from the photometric age- 
dating technique that results in an underestimation of 
the number of clusters in this age interval. It is com- 
monly claimed as being an artifact from the photometric 
age-d ating technique (e.g. iGieles et al.l 120051 iLee et al.l 
I2005f ) , but more recently iMaiz Apellanizl (|2009T ) pro- 
posed that it is more specifically related to stochastic 
effects in the stellar initial mass function sampling. In 
Figure [5] we also indicate for each band the magnitude 
of the brightest stars, and where confusion can occur be- 
tween a bright star a nd a faint cluster. T he limits are 
based on the work of iMassev et al.l (|2006l their Fig. 4 
to 6) that includes a large sample of stars in M31 and 
M 33. In the V ban d, this is consiste nt with the work 
of lHumphrevsl (|1983f l. Iffill et al.l (|1993l ) obtained a good 
stellar photometry in the B-band for massive stars in 
30 Do r and reporte d a lim it value that is slightly fainter 
than IMassev et al.l (|2006D . Therefore an averaged limit 
value was used in this work for the B band. Note that 
clusters more massive than lxl0 5 M Q do not fall below 
the star/cluster brightness limit for UBV bands. 

Age uncertainties are displayed in Figure |6] as a func- 
tion of cluster brightness and age. It shows no partic- 
ular trend other than larger uncertainties are observed 
for fainter sources. Square symbols indicates sources for 
which the completeness level is lower than 50%. Note 
that the BVI sample includes only such sources. The 
bottom panel clearly shows the population of old clusters 
(>300Myr) has a low completeness level, and is therefore 
not reliable to draw scientific conclusions. 

To determine the age at which the sample is com- 
plete, we use Figure [5] and find the age where the model 
lines (solid) cross the lower envelope of the data points 
with >50% completeness and the horizontal (dashed) line 
marking the brightest stars. The lowest age of these two 
determines the age completeness for that filter and mass 
limit, and is marked by a dotted line in each panel of 
Figure the lowest of these for U, B, V, and I deter- 
mines our final age limit, since clusters must be detected 
in at least these filters to have a photometric age. This 
results in an age limit of 125 Myr for the 1 x 10 5 M Q sam- 
ple which is effectively defined by both the brightest star 
luminosity and the 1 x 10 5 M model line in the U band. 
Those values indicate that we do not have to worry about 
statistical fluctuations of the initial mass function that 
can be s een for clusters less massi ve than a few times 
1O 4 M (ICervino fc Luridianal2004h . Also, it means that 
we cannot extend our mass-limited sample analysis to 
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Fig. 3. — Cluster luminosity function for NGC922. Black: all sources detected in each individual filter. Blue: Sources detected 
simultaneously in the UBV1 filters. Red: Sources detected simultaneously in the BVI filters only. Green: combination of the UBVI plus 
BVI samples (i.e. Blue+Red curves). The dotted lines represent the samples uncorrected for incompleteness. The arrows show the 50% 
completeness limit. 



ages greater than ~ 125 Myr. 

The age distribution of stellar clusters in NGC 922 is 
presented in Figure [7] A detailed list of the completeness 
corrections applied to Fig. [7J is presented in Table El The 
age interval from 1 to 330 Myr covers essentially the time 
over which NGC 922 should have experienced a relatively 
high star formation activity and produced a large amount 
of star clusters due to the collision with S2, and the sub- 
sequent propagation of the density wave through the disk 
(W06). The number of clusters was corrected for incom- 
pleteness (see $2]), an d the age distributions include only 
sources with combined completeness levels higher than 
50%. 

We find slopes of — 0.7±0.2 for the age distributions 
in Figure [7J The slopes were obtain ed using ordinary 
least- squares of Y|X linear regressions (|Feigelson fc Babul 
[1991) over the mass ranges that are not affected by in- 
completeness, i.e. 6.5 < log(r) < 8.0 if using a bin size 
of 0.5 dex and 6.6 < log(r) <8.1 if using a bin size of 
0.3 dex. Note that a bin size of 0.3 dex is smaller than 
the typical age uncertainties as shown in Fig [5] We ex- 
cluded the first bins (0-5 Myr) because clusters may still 
be buried in their native dust cocoon at this young age, 
adding to the incompleteness at visible wavelengths. We 
also excluded the bin around 30 Myr in Figure [JJd due 
to the age bias discussed above. The 30 Myr bin was 



included in the fit in Figure [JJi because of the very low 
number of bins available for the fit. We also performed 
a Kolmogorov-Smirnov (KS) test as an independent test 
to check the slope of the observed age distribution. The 
test is presented in Fig. [8] and show results consistent 
with negative slopes. We find a probability of less than 
1% for a positive slope, while we obtain the highest prob- 
abilities (>50%) for slopes between —0.4 and —0.8. Note 
that due to jumps in the cumulative distribution, we can- 
not find a single best fit. 

The evolution of the SFR as traced by star clusters 
should, to first order, tell us something about the distri- 
bution of the dense interstellar medium (ISM), or more 
specifically the molecular gas, from which the clusters 
formed. In a collisional ring galaxy much of the star for- 
mation occurs in the expanding ring. Hence, one may 
expect the SFR to be proportional to the rate of the 
dense ISM being swept up by this ring. If the expan- 
sion velocity is nearly constant, and the ring radius is 
significantly larger than the radius at which the impact 
occurred, and if there is no cluster dissolution then a 
flat cluster age distribution would mean that the dense 
ISM must have a surface density S g that falls off with 
radius, r, as r _1 . Wh ile the neutral ISM in galaxies can 
fall off this shallowly (|Bosmalll98illHoekstra et al.|[200ll: 
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Fig. 4. — Maps of star clusters detected in NGC 922 for both UBVI and BVI samples, a) Clusters younger than 7Myr and capable of 
producing Ha emission, b) Young clusters with ages from 7 to 30Myr. c) Clusters with ages from 30 to 150 Myr. d) Clusters older than 
150 Myr. The grey star symbols show the location of the two unresolved clusters and the nucleus (see Table [2j|. 



TABLE 3 

Incompleteness Corrections of Observed Cluster 
Numbers and Age Distribution Values 



Age Bin 


No. Clusters 51 


No. Cluster 13 


Log(dN/dt) 




Uncorrected 


Corrected 


Corrected 


0.3 dex bins 








6.6-6.9 


19.0±4.3 


23.7±5.4 


-5.22 


6.9-7.2 


8.0±2.8 


9.8±3.5 


-5.91 


7.2-7.5 


4.0±2.0 


4.7±2.3 


-6.53 


7.5-7.8 


26.0±5.1 


30.9±6.1 


-6.01 


7.8-8.1 


32.0±5.7 


39.5±7.0 


-6.20 


0.5 dex bins 








6.5-7.0 


20.0±4.5 


24.7±5.5 


-5.44 


7.0-7.5 


11.0±3.3 


13.5±4.1 


-6.21 


7.5-8.0 


43.0±6.5 


52.3±8.0 


-6.12 



a Number of clusters observed and uncorrected for incomplete- 
ness. The values corresponds to the dotted lines in Fig. fTb.d. 
k Number of clusters observed and corrected for incompleteness. 
The values corresponds to the solid lines in Fig. [7b, d. 
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Log ir/yr) Log (r/yr) 

Fig. 5. — Magnitudc-agc distributions of star clusters observed in NGC922, including both UBVI and BVI samples shown in term of 
their a) U, b) B, c) V, and d) I magnitudes. Solid circles: clusters with combined completeness level greater or equal to 50%. Open 
squares: clusters with combined completeness level smaller than 50%. The solid lines represent stellar populations with masses of lxlO 4 , 
1x10 s , and lxlO 6 Mq, from bottom to top. The vertical dotted lines indicate the age at which the sample of clusters more massive than 
lxlO 5 M0 is complete in the specified filter. The horizontal dashed lines determine the brightness upper limit for a single star in each 
band. The right axis is displayed in absolute magnitude. 
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Fig. 6. — Top panel: Age uncertainties as a function of cluster brightness. Bottom panel: age uncertainties as a function of cluster age. 
Solid circles: clusters with completeness level greater or equal to 50%. Open squares: clusters with completeness level smaller than 50%. 
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Fig. 7. — a) Observed age distribution of star clusters of at least 1x10 s Mq using a bin size of 0.5 dex in log(age ). b) Same as in panel 
a) but for a bin size of 0.3 dex. Linear fits for 6.5>Log("r/yr)>8.0, where the sample is considered as complete (see !|4,2jl . are over-plotted, 
c) Simple cluster age histograms for the same sample using a bin size of 0.5 dex. d) Same as in panel c) but using a bin size of 0.3 dex. 
Dotted lines: sample uncorrected for incompleteness. Solid lines: sample corrected for incompleteness. Clusters with completeness lower 
than 50% are excluded. 



iLeroy et al.ll2008T ). the molecular ISM, from which stars 
form, has a steeper fall-off m ore similar to the stellar 
distribution (|Lerov et al.l 120081) . This will produce posi- 
tive slopes in log(iV) versus log(age), as we indeed find 
in the simulations presented in Section 6. A negative 
slope suggests a shallower drop-off than £ g (V) oc r _1 ; 
while one can not rule this out in any particular galaxy, 
it would represent a peculiar molecular ISM distribution. 
Assuming a normal steep S s (r) profile, one way to rec- 
oncile the negative slope observed in NGC 922 with what 
we understand of collisional ring galaxies is to introduce 
young star cluster disruption and dissolution. In the next 
section we explore this possibility. 

5. STAR CLUSTER DISRUPTION IN NGC 922 

The disruption of star clusters is a concept not triv- 
ial to test observationally, especially because the star 
for mation history is oft en unknown. The recent work 
of iBastian et all (|2009h expresses well this problem. 
NGC 922 represents a good case-study for the star clus- 
ter disruption scenario since its cluster formation history 
can be relatively well recovered with dynamical and star 
formation simulations. 

In t heir work on the Antennae galaxies, [Fall et aT] 
( 2005) proposed that about 90% of star clusters dissolve 
for each dex of age, up to 1 Gyr, independently of clus- 
ter mass. We will call this empirically based relation 
the Mass Inde pende nt ( MI) scenario, following th e works 
of iFall et all ([2005^ and iWhitmore et al.l ([200l . Con- 
cerns about the mass independent, and continuous na- 
ture of the cluster destruction, as well as t he assumption 
of a constant SFR have been raised (e.g. IBastian et al.l 
2009). There is an alternative scenario in which clus- 
ter dissolution rate varies over time. According to this 



scenario, infant mortality is responsible for cluster dis- 
solution during the first 10 7 yr or so, leading to a steep 
decreasing slope in the age-distribution. Infant mortal- 
ity of star clusters is a phenomenon thought to hap- 
pen to young clusters experiencing a serious loss of gas 
and dust during the supernovae explosion phase, leav- 
ing the clusters gravitationally un bound and easy to 
disrupt within 10 to 50 Myr ( e.g. ILada fc Ladal 120031: 
IFall et all [20051 ; iGoodwinl 12001) . According to a recent 
work, i nfant mortality is hard ly seen at visible wave- 
lengths (Gieles & Bastian 2008). Then internal dynam- 
ical evolution and stellar population fading occurs be- 
tween 10 7 and 10 8 yr. There are relatively few clus- 
ter destroyed during that phase. This results in a flat 
age distribution as sketched by Lamers (2008), although 
the precise form may not be exactly flat, which he also 
notes. Finally mass dependent cluster dissolution, due 
to galactic tidal effects and other external effects, pro- 
duces a decreasing slope for the 10 8 -10 9 yr. We will call 
this scenario the Mass Dependent ( MP) sc enario, after 
the w or ks of iBoutloukos fc Lameri (I2003T), IGieles et al.l 
( 2005) . Il^rners et ; all JI2005D . ILamers fe Gielesl (|2006) , 
and IGieles et all (|2007l) . A good compari s on be tween 
the two scenarios is presented by ILamers! ((2008). As 
shown in Figure [3 there might be some evidence of a 
flat distribution in NGC 922 between 10 and 100 Myr, as 
proposed by the MD scenario, but that could also be due 
to the artificial lack of clusters around 10-30 Myr men- 
tioned above. 

For every cluster detected in NGC 922, we applied a 
correction for the loss of stellar clusters due to dissolu- 
tion. We used the two cluster disruption functions dis- 
cussed above. A starting age of 4 Myr was chosen in- 
stead of Myr because it corresponds better to the age 
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Fig. 8. — Cumulative distribution function of star clusters of at least 1x10 s Mq (circles). For comparison, the dashed, dotted, plain, and 
dash-dotted lines show what would be the cumulative distribution function of an age distribution with slopes of — 1, —0.6, 0, and +0.37, 
respectively. The positive slope represents a typical case of a ring galaxy the size of NGC 922 as simulated in 36.31 



of the first supernova explosions and, presumably, the 
time when the clusters can start dissolving. The fraction 
of star clusters F sc left after dissolution, according to the 
MI model, is then given by 

Zo 5 (F sc ) = 1, for log(r) < 6.6 (1) 

log(F sc ) = 6.6 - log(r), for log(r) > 6.6 (2) 

where r is the age of the cluster in years. Following the 
MD model, the fraction can be expressed by 

log(F sc ) = 1, for log(r) < 6.6 (3) 

log(F sc ) = 16.5 - 2.51og(r), for 6.6 < log(r) < 7.0 (4) 

log(F sc ) = -1, for 7.0 < log(r) < 8.0 (5) 

Zo 5 (F sc ) = 12.6 - 1.71og(r), for log(r) > 8.0 (6) 

For the early ages of the MD model, we used a 90% 
infant mortality rate per decade of time. Note that 
because we do not know the star formation history of 
NGC 922 prior to the collision, we cannot discuss the is- 
sue of whether or not the cluster disappearance in this 
age range is mass-dependent (see i )4.2[) . 

The age distributions of star clusters corrected for clus- 
ter dissolution and incompleteness is presented in Fig- 
ure IH1 We find positive slopes of +0.5±0.2 for clusters 
with masses of at least 1 x 10 5 M Q when we use the MI 
model. We find a relatively fiat slopes (0.0±0.1) when 
we use the MD model. The slopes were calculated up to 
Io 9 {t) = 8.1. 

The positive slope value below 100 Myr now sug- 
gests that the cluster formation rate was decreasing in 
NGC 922 during that period, assuming that the dis- 
solution correction is realistic and that there are no 
other major mechanisms affecting the number of clus- 
ters. The slope value for clusters older than 100 Myr is 



still steep despite the correction for dissolution but this 
is attributed to incompleteness. 

6. DYNAMICAL AND STAR-CLUSTER SIMULATIONS 

In order to better understand the age distribution of 
star clusters observed in NGC 922, we performed N- 
body/SPH simulations on possible formation sites of star 
clusters in NGC 922 colliding with its companion S2. 
We considered that the observed ring-like morphology 
of NGC 922 is the result of a high-speed, off-centered 
collision between NGC 922 and S2 (W06). We therefore 
describe the simulation results with unique orbital con- 
figurations and disk inclinations. Here we focus mainly 
on the age and spatial distributions of star clusters within 
NGC 922 and not the chemical and morphological evolu- 
tion. 

6.1. Model Parameters 

In the current simulation NGC 922 is represented by a 
fully self-consistent disk galaxy model whereas S2 is rep- 
resented by a point mass. NGC 922 is modeled as a disk 
galaxy embedded in a massive dark matter halo with 
a total mass of Mdm and the universal "NFW profile" 
([Navarro et al.lll996r ). The disk has an exponential radius 
and the total disk mass and size are M<j and Rd, respec- 
tively. Our results were obtained using =2x 10 10 M Q , 
i?d = 13.4 kpc, and M^m/Md — 9. The gas mass fraction 
/ g of the disk is assumed to be a free parameter. More 
details about the numerical methods and techniques we 
employ for modeling chemodynamical evolution of in- 
teracting ga l axies with star formation can by found in 
IBekki et all (pOOl . 

The radial (R) and vertical (Z) density profiles of the 
disk are assumed to be proportional to exp(— R/ Rq) with 

scale length Rq = 0.2i?d and to sech 2 (Z/Zo) with scale 
length Z — 0.04i?d, respectively. In addition to the 
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Fig. 9. — a) Observed age distribution of star clusters of at least lxlO 5 M0 and corrected for cluster disruption and incompleteness, 
b) Same as in panel a) but for a 0.3 dex bin size, c) Simple cluster age histogram for the sample corrected for cluster disruption and 
incompleteness, d) Same as in panel c) but for a different bin size. The solid lines show the sample corrected using the MI model while 
the dashed lines show the sample corrected using the MD model. The error bars on the MD corrected histograms are omitted here so that 
the plots are not too cluttered. They are similar in size to those for the MI corrected histograms. 



rotational velocity attributable to the gravitational field 
of the disk and halo components, the initial radial and 
azimuthal velocity dispersions are added to the disk com- 
ponent in accordance with the epicy clic theory, and with 
a Too mre Q parameter value of 1.5 (|Binnev fe Tremaind 
119871 ). The interstellar medium in NGC 922 is repre- 
sented by SPH particles. The mass resolution (i.e., the 
minimum mass of a particle) and the size resolution (the 
minimum gravitational softening length) are 2xl0 5 M© 
and 13.4 pc, respectively. 

The mass ratio of the companion dwarf to the spiral 
(7712) is set to be 0.2. The initial location and velocity 
of the companion S2, with respect to be the center of 
the spiral, are represented by X g and V g , respectively. 
The spiral is assumed to be inclined by 80 degrees with 
respect to the x-y plane. Although we have investigated 
models with different X g and V g , we show the results of 
the model with (x,y,z) = (— 4Pd, 0.2Pm, 0) and (v x ,v y ,v z ) 
= (6v, 0, 0), where v is the simulation unit (= 80.1 km 
s^ 1 ). Stars that are initially within the spiral disk are 
referred to as "stars" (or "old stars" ) whereas stars that 
are newly formed from gas are as "new stars" . 

6.2. Possible formation sites of star clusters 

We adopted the following star formation laws to con- 
vert gaseous particles into either star clusters or field 
stars. For star cluster formation, we considered that in- 
terstellar gaseous pressure (P ga s) in star forming regions 
of a galaxy drives the collapse of pressure-confined mag- 
netized self-gravitating molecular clouds to form compact 
clusters, providing that P gas is l arger than the surface 
press ure (P B ) of the clouds (e.g. lElmegreen fc Efremovl 
119971) : 



Pgas > P s ~ 2.0 x 10 5 fc B . (7) 

where ks is the Boltzman constant. This adopted value 
for the threshold gaseous pressure P gas would be appro- 
priate for the formation of massive star clusters, such as 
the lxl0 5 M Q mass-limited sample. Our star formation 
prescription allows for the formation of low mass clusters, 
associations or field stars even if P gas is smaller than the 
above threshold value. However, we do not track these 
stars in the plots presented here. 

For field stars, we adopt the Schmidt law dSchmidtJ 
119591 ) with an exponent of 1.5 (jKennicuttl Il989f ). Al- 
though we cannot investigate the detailed physical pro- 
cesses of cluster formation in the present and global sim- 
ulation (i.e. from ~ 100 pc to 10 kpc scale), we expect 
that the adopted phenomenological approach enables us 
to identify the plausible formation sites of star clusters 
in NGC 922. For models with / g = 0.4 (note that this 
is different from the 0.2 value adopted by W06), the 
present possible mass of young stars in NGC 922 can be 
well reproduced. For the adopted orbital configurations, 
NGC 922 can experience a moderately strong starburst 
when S2 passes through NGC 922. 

6.3. Simulation Results 

The best-fit model produces a peak of star formation 
during the passage of S2 within the disk of NGC 922 
about 150 Myr ago, then a relatively constant decrease of 
the cluster formation rate over time. Figure [TUl shows the 
simulated age distribution of 4240 star clusters formed 
as described above (see a detailed discussion about the 
star cluster or g lobular cluster formation models by 
iBekki et al.l 120021) . The overall shape of the simulated 
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TABLE 4 
Corrected 2MASS Photometry for 
NGC922 and S2 



NGC 922 



S2 



NUV (AB mag) 
R (AB mag) 
J (ABmag) 
H (AB mag) 
K (AB mag) 
M* (M Q ) 



13.04 
11.57 
11.57 
11.49 
11.82 



15.98 
14.71 
14.87 
14.89 
14.80 



2.8xl0 10 1.81x10 s 



Note. — All values have been corrected 
for extinction as in W06. 

age distribution displays a slope of +0.37±0.02 between 
1 and 100 Myr. 

The expected cumulative distribution function of clus- 
ters formed with this age distribution and no disruption 
is shown in Fig. |8] with a dot-dashed line. It shows that 
the age slope of —0.7 ± 0.2 we find without any clus- 
ter dissolution correction is significantly steeper than our 
simulation result by ~ 5a. As discussed in Section [5] and 
illustrated in Fig. |H1 after correction for MI we find an 
age slope of 0.5 ± 0.2 which is within ~ 1 a of the simu- 
lated results. Correction using the MD scenario produces 
a flatter age slope of 0.0±0.1, which differs from our sim- 
ulation by ~ 4cr. We therefore conclude that correction 
for cluster dissolution, using previously defined prescrip- 
tions, is capable of bringing the observed age distribution 
into better accord with the simulations, and that the MI 
model for the correction produces results somewhat more 
consistent with the simulations than the MD based cor- 
rection. Because of the possibility that any given galaxy 
may have a peculiar cluster formation history that is not 
captured in a simple model, one should not see our re- 
sults as a sweeping proof of cluster destruction following 
the MI scenario, the MD scenario still agrees within 4cr 
of the simple model. However, it does suggest that fur- 
ther comparisons between simulations and well modelled 
galaxies may provide good constraints on which cluster 
disruption path works best. 

One possible problem with the MI scenario is the vast 
number of clusters that would have had to be destroyed. 
To test if the implied SFR is consistent with the obser- 
vations of W06, we integrated the MI cluster formation 
rate for all detected clusters uncorrected for incomplete- 
ness. If we use the MI scenario, we obtain a total stellar 
mass post-collision of 7x 10 8 M Q . This is 14% of the total 
stellar mass reported by W06 (5x10 9 Mq) for the whole 
galaxy. This stellar mass is a lower limit if we consider 
that the collision was 300 Myr ago as suggested by W06, 
and that our sample is not complete beyond 125 Myr 
ago. However, we found an error in the conversion of 
the 2MASS fluxes to ABmag presented by W06. We re- 
calculated the photometry and masses (see Table [4]) and 
obtain a total stellar mass of 2.8xlO lo M for NGC 922 
based on the K-band, which makes the corrected mass 
originally created in clusters to be 2.5% of the total stel- 
lar mass. We concluded that the integrated cluster mass 
produced in this starburst is not an unreasonable fraction 
of the total mass. 

The simulated star clusters younger than 100 Myr are 
located preferentially in the ring. This result confirms 
that the high-density, high-pressure ring formed after 



the NGC 922-S2 collision is the ideal site for star clus- 
ter formation, and is consistent with the observations. 
Also, according to the present simulation, NGC 922 and 
the companion collided about 150 Myr ago which trig- 
gered a starburst in the central region of NGC 922. This 
burst population does not show a ring-like distribution 
and is thus responsible for a more dispersed distribu- 
tion for star clusters with ages between 100 and 200 
Myr. This population is not obvious in the HST im- 
ages, but this is related to the detection limit of our 
sample or cluster dissolution, most of the older clus- 
ters being seen either within the nuclear region or the 
ring. However, in contradiction to the models, the 
observations show that star clusters also formed in- 
side the ring within the last 100 Myr. Those were 
possibly formed along some spoke f eatures, as shown 
in other collisional ring simulations ([Hernauist fc Weil 
fll^lAmHeton fc Struck-Maroel]llT99fiHStrucldll997i r~ 

7. LARGE-SCALE STAR FORMATION 

Large-scale star formation in galaxies provides in- 
sight into global dynamical processes. In normal spi- 
ral galaxies, the largest star-forming units are kpc-size 
complexes which incorporate many smaller clumps or 
clusters. The sizes of these largest complexes scale 
with th e absolute magnitude and the diameter of the 
galaxy ([Elmegreen et all 11994 he reafter EELS). Th ey 
are governed by a Jeans instability (jElme grecn 1994bllal). 
which oft en leads to rather regular sp acing along spi- 
ral arms (Elmegreen & Elmegreen 1983) and in circum- 
nuclear rings (iButa fc Crocker! 119931: iMaoz etaLl 119961: 
lElmegreen et al.l I1997D . The gas mass of H i super- 
cloud complexes is on the order of 10 7 M Q for spiral 
arms (iElmcgrc en fc Elmegreenlll987D and the total stel- 
lar mass upper limit is less than 10 6 M consisting of 
several embedded clusters; in circumnuclear rings, the 
hotspot stellar masses are typically a few times 10 5 M Q 



(|Elmegreen et al.lll997j ) . Complexes are also the smallest 
structures that can be resolved in high redshift galaxies, 
so these star-forming units are useful to examine in NGC 
922 for comparison with ring galaxies in the early uni- 
verse. 

7.1. Star-forming complexes 

Complexes in NGC 922 were identified by visual in- 
spection of the WFPC2 images. Figure [TT] shows the 
F547M image with the star-forming complexes labeled; 
these encompass the smaller-scale clusters discussed in 
sections 2] to IH1 These complexes, which range from 750 
pc to 2 kpc in diameter, are also prominent in the SINGG 
Ha image of NGC 922 (see Fig. [TJ. In nearby quiescent 
spiral galaxies, a late-type gal axy with NGC 922' s abso- 
lute magnitude of M fl =-21.24 ([Meurer et al.ll2006[ ) would 
be expected to have a maximum complex size of ~1 kpc 
(see EELS). 

Photometry was done on the complexes in each of the 
WFPC2 images by defining boxes around the complexes 
and their brightest subunits and using the IRAF task 
imstat. The results were compared with evolutionary 
models of an instantaneous burst, as in section 3. The 
stellar population synthesis results are presented in Ta- 
ble [5l The complexes have ages ranging from 4 to 86 Myr 
and masses from 9xl0 5 to 3xlO 7 M . 
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Fig. 10. — Top: simulated age distribution of clusters formed after NGC 922's collision with S2. Bottom: Simulated cluster age histogram. 



TABLE 5 

Physical Parameters of Star-Forming complexes 



Complex ID 




RA 


DEC 


F300W 


F439W 


F547M 


F814W 


Age 


Log(M Ste ((ar) 


Diameter 




degrees 


degrees 


mag 


mag 


mag 


mag 


Myr 


M 


kpc 


1 


36 


272875 


-24.798485 


15.49±0.07 


16.05±0.07 


16.6±0.1 


17.4±0.2 


4 +B0 


6 7+ 1 ' 1 
D - ' -0.3 


0.76 


2 


36 


275468 


-24.795138 


16.72±0.04 


17.01±0.03 


17.46±0.08 


18.1±0.1 


7+92 


D.O_ Q 2 


0.71 


3 


36 


275073 


-24.790066 


16.00±0.01 


15.97±0.01 


16.42±0.03 


17.02±0.04 




7.6+J- 5 


1.04 


5 


36 


274900 


-24.786761 


17.39±0.02 


17.67±0.02 


18.13±0.05 


18.70±0.07 


7+291 


6 0+ 2 ' 3 
°- u -o.i 


0.43 


6 


36 


274171 


-24.783779 


17.17±0.01 


17.25±0.02 


17.65±0.04 


18.19±0.07 


c 4 +303 


7.1+1' 4 


0.6 


8-9 


36 


272113 


-24.784431 


16.17±0.02 


16.13±0.02 


16.52±0.04 


17.05±0.07 


54+50 




0.92 


11 


36 


266513 


-24.800236 


16.63±0.02 


16.96±0.02 


17.45±0.04 


18.07±0.06 


7+8O 


6.3±S"i 


0.78 


12-13 


36 


268845 


-24.801279 


16.78±0.02 


17.02±0.02 


17.52±0.04 


18.13±0.06 


71?™ 


6.212-3 


0.76 


14 


36 


270572 


-24.800231 


16.77±0.01 


17.02±0.01 


17.53±0.03 


18.14±0.05 


7+290 


6 2+ 2 ' 4 
D - z -0.3 


0.71 


15 


36 


256402 


-24.797122 


16.34±0.01 


16.47±0.01 


16.85±0.03 


17.43±0.04 


54+268 


7 4+1.4 
'■^-1.2 


1.19 


16 


36 


264253 


-24.794724 


16.26±0.005 


16.05±0.008 


16.47±0.02 


17.05±0.02 


54+«° 


7 5+ 1 ' 5 
'■°-1.3 


1.19 


17 


36 


265672 


-24.793388 


16.08±0.005 


16.52±0.008 


16.88±0.02 


17.42±0.03 


87+4? 


7 5+ 1 ' 4 
'■°-1.6 


0.91 


18 


36 


262733 


-24.790521 


17.51±0.005 


17.36±0.007 


17.72±0.02 


18.28±0.02 


54t|37 


7Q+1.6 
'■"-1.4 


0.67 



Note. — The RA and DEC positions arc in J2000 and approximations and arc not necessarily related to the exact photometric center of the 
sources. 
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7.2. Comparison with High Redshift Collisional Ring 

Clumps 

The largest complexes in NGC 922 are close to the limit 
of resolved star-forming re gions in galaxies in th e Hub- 
ble Ultra Deep Field (UDF iBeckwith et al.|[2006h out to 
redshifts of about 4-5; W06 n oted that several clump 
clust er galaxies in the UDF (jElmegreen fe Elmegrcen 
|2005[ ) resemble NGC 922 in having outer rings. Galax- 
ies in the GEMS and GOODS fields at redshifts 0.6 
to 1.4 also include possible collisional ring galaxies 
(jElmegreen fc E lmegreenl I2006D . Two UDF galaxies, 
with catalog numbers 5190 and 9159, have large-scale 
star-form ation in their outer ring-like structures (see 
Fig. 1 in lElmegreen fc Elmegreenl 120051 ). Their photo- 
metric redshifts are 1.25 and 1.45, respectively. They 
each have a diameter of about 10 kpc, and the com- 
plexes range from 0.5 to 2 kpc. The largest complex 
is about 13% of the galaxy size in UDF5190 and 20% 
in UDF9159, which is 3 times larger than the average 
for nearby late-type spirals but about the same as the 
ratio for NGC922, 15% based on its diameter of 13.4 
kpc (W06). The absolute complex sizes are about the 
same in the distant galaxies as in NGC 922, but the dis- 
tant galaxies are physically smaller. Photometric data 
for the clumps in these two galaxies were compared with 
evolutionary models that included bandshifting effects, 
interstellar re ddening, and intergalactic Hi abso rption 
(see details in lElmegreen fc Elmegreenll2005L |200H . The 
clump ages are about 10 8 yr, and their masses range from 
2xl0 7 to 2xl0 8 M Q . Thus, they are fairly similar to the 
largest complexes in NGC 922. 

A GOODS galaxy with a well-defined outer ring of star 
formation and an off-center nucleus is shown in Figure 
[T2l (left) in color from the MAST archivetQ. It is number 
21238 in the COM BO-17 catalog, with a redshift of 0.664 
(Wolf ct al. 2003). Star- forming complexes are identified 
in Figure H"2l (right). Number 1 is the nucleus; 2 is the 
brightest star-forming region on the ring, and other ring 
clumps are 3-13. Features 14-16 are "feathers" on the 
western side, reminiscent of the plume in NGC 922, and 
17 is a blue clump. Photometry was done on the fea- 
tures in the BVIz images. Figure [13] shows their color- 
magnitude diagram of V606 versus (V606" z 85o)- Over- 
laid are evolutionary tracks for clusters with different 
masses, as labeled. The open circles along the lines are 
timesteps, with two times indicated. Most of the clumps 
are 3-7xl0 7 M Q . The feathers are in the same range. 
The blue dot, Feature 17, is much less massive, closer to 
lxlO 7 M . The nucleus is about 5 x10 s M Q . The ages 
of the ring clumps are 80 to 150 Myr. The clumps on the 
southeast side of the ring are systematically redder than 
those on the northwest, possibly due to some reddening. 

Other GEMS and GOODS ring galaxies at redshifts 
z=0.7 to 1.4 have ring clumps that are about 0.1 
times the galaxy diameter, with ages of a fewxl0 8 yr. 
Their masses range from several x 10 6 M Q for the 
lower redshifts to a fewxlO 8 M^for the higher red- 
shifts (|Elmegreen fc Elmegrecn 2006|). These masses 
are slightly smaller than the giant clumps found 
in several non-ring clump cluster galaxies at red- 
shift z=1.8 to 3.0, whose masses average 6x10 s M Q 

6 http: /archive. stsci.edu 



(Elmegrcen & Elmegrcen 2005). 

8. THE TIDAL PLUME AND DIFFUSE EMISSION OF 

NGC 922 

The three-color image of NGC 922 (Fig. [J) shows a 
faint yellow-red tail on the western side of the nucleus, 
pointing roughly toward the companion S2. This struc- 
ture appears clearly in the WF2 chip. A smooth and 
unresolved stellar population comprises the tail; any clus- 
ters present are below the detection limit in U and B. 

As was done for the nucleus and the two large clusters 
and nucleus (Section [2|), we made photometric measure- 
ments for the tail (see Table [2]) in the area shown by 
a polygon in Figure [2] There was no detection in the 
F300W image. The photometric uncertainties depend 
mainly on the size of the polygon. 

Our population synthesis of the plume indicates stellar 
ages of 1 to 5 Gyr, and a total stellar mass of the order of 
7xl0 8 M Q . The age found for the tail strongly suggests 
either that it was part of the main disk of NGC 922, and 
was pulled away by S2 after the collision or it was part of 
S2 and was ripped off during the encounter. Based on the 
I-band residual image, we can see the tail extending up to 
7-8 kpc from star forming regions in the disk. If the stars 
come from NGC 922's disk, this implies that the most 
distant stars have a projected linear velocity of about 
20-30 km sec -1 . Even though the material in the tail was 
subject to intense mixing and compressions, it did not 
produce significant star formation. This suggests that 
the plume may have formed from an ISM-frcc portion of 
the galaxy (i.e. bulge or thick disk). 

In addition to the tidal plume, another diffuse emis- 
sion component is observed in the whole disk of NGC 922 
after removing the clusters from the images. We calcu- 
lated the fraction of diffuse light and found that 79.2, 
89.5, 93.0, and 94.1% of the light was produced by this 
diffuse component for the U, B, V, and I filters, respec- 
tively. This indicates the fraction of light in clusters is 
more important at shorter wavelengths. The fraction 
values are also consistent with the work of iMeurer et al.l 
(1995 !) where they found that about 80% of the UV light 
is found outside of star clusters in nearby starburst galax- 
ies. The diffuse fractio n is known to be higher in normal 
galaxies (|Larsenll200Tjl ). This indicates that NGC 922 is 
more a starburst type of galaxy, i.e. with enhanced clus- 
ter formation activity. Synthesis modeling suggests an 
age of 250±40Myr for the diffuse emission, indicating 
that the diffuse stellar component was mainly formed 
before the collision ^150 Myr ago. We obtained that 
log(Ms te i2 ar ) = 10.0 ± 0.1 M© for the diffuse emission, 
although the use of an instantaneous burst model might 
be less appropriate in this case. 

9. THE FATE OF NGC 922 AND ITS STAR CLUSTERS 

According to general dynamical models, the ring in 
NGC 922 will eventually be fragmented into several mas- 
sive clumps, which would resemble high-redshift com- 
plexes and substructures. Its galactic disk will become 
thicker and it will evolve in morphology to become an 
SO. 

On a smaller scale point-of-view, iBurkert et all (|200l 
proposed that collisional ring galaxies were excel- 
lent candidates for producing faint fuzzy clusters. 
They are known to be old, relatively large and non- 
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Fig. 12.— Left: Collisional ring galaxy COMBO-17 21238 from the Hubble GOODS field number 12. Right: The complexes and features 
are labeled. See section [7.21 for a detailed identification of the morphological features. The image is 3" in length, or about 25kpc. 
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Fig. 13. — Color-magnitude diagram showing the complexes and feathers in COMBO-17 21238, with evolutionary tracks overlaid showing 
cluster ages and masses. 



compact clusters (8-15 Gyr , 7-15 pc lLarsen fc Brodid 
[200aiBrodie fc Larserj|2002l) . lEhnegreenl (|2008D showed 
that they are actually the lucky survivors of star forma- 
tion in regions where the gas density was barely enough 
to produce bound clusters but where the environment 
also presented low tidal forces so the clusters were not 
torn apart. The author explains that such environments 
can exist in the molecular clouds of outer galactic regions 
that are highly shocked such as in the expanding colli- 
sional rings of star formation. Consequently, the most 



massive older clusters (e.g. clusters >50Myr in Fig. [5]) 
observed in NGC 922 are excellent candidates for progen- 
itors of faint fuzzy clusters. 

10. CONCLUSIONS 

In this work, we presented the four-band HST/WFPC2 
images of NGC 922, one of the nearest known collisional 
ring galaxies. We performed PSF photometry and de- 
tected 3817 stellar clusters in at least three bands, as 
well as two other clusters bigger than the PSF. We used 
the stellar population synthesis model Starburst99 to es- 
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timate ages and masses for each observed clusters. We 
also simulated the dynamical and star-forming events re- 
lated to the collision between NGC 922 and its compan- 
ion S2. 

We found that most observed clusters are 7Myr old, 
or younger. Those clusters are mainly located within 
the ring or the bar, consistent with Ha emission and in 
agreement with the simulation. A large fraction of clus- 
ters older than 50Myr are found in the nuclear region, 
which is consistent with the model producing many older 
clusters (100-200 Myr) in the bulge. Those also tend to 
be massive, according to single stellar population model- 
ing, but are probably contaminated by an older underly- 
ing stellar population which would then bias the age and 
mass estimation. 

We also observed a stellar tidal plume pointing in the 
direction of the companion. We did not detect clusters in 
the plume. We performed integrated photometry for this 
plume and found that it is consistent with stellar popu- 
lations of 1 to 5 Gyr old. This suggests that the material 
was stripped off from NGC 922 by S2 and did not pro- 
duce significant star formation related to the collision. 

The cluster luminosity function observed for NGC 922 
has a slope typical of young cluster populations. The 
age distribution of star clusters shows a decreasing slope, 
consistent with cluster dissolution and disruption. By 
comparing the N-body/SPH star formation simulation 
to the observed age distribution corrected for cluster dis- 
ruption (MI and MD models), we find that the MI model 
is over-correcting for dissolution while the MD model is 
under-correcting in the 10-100 Myr age range. In any 



case, cluster dissolution and disruption is a plausible 
physical process to consider in order to reach good agree- 
ment between the simulation and observation. 

We found a significant number of massive (>1O 5 M ) 
star clusters in NGC 922 that are old enough (>50Myr) 
to have survived the supernova explosion phase. Since 
some of them are likely to be born in highly shocked 
clouds and that they are now located at high galactic 
radii where the tidal disruption forces are presumably 
low, those clusters are excellent candidate progenitors of 
faint fuzzy clusters. 

NGC 922 has a morphology similar to several distant 
galaxies observed in the GEMS and GOODS fields. We 
detected about a dozen big star-forming complexes in the 
ring of NGC 922. A direct comparison with the knots ob- 
served in a ring galaxy from the GOODS field (#21238 in 
COMBO-17, z=0.7) and other high redshift ring galax- 
ies shows that NGC922's complexes are of comparable 
sizes and masses even though the distant galaxies are 
physically smaller. 
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